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We investigate a condensed matter "black hole" analog, taking the Gross-Pitaevskii (GP) equation 
as a starting point. The linearized GP equation corresponds to a wave equation on a black hole 
background, giving quasinormal modes under some appropriate conditions. We suggest that we can 
know the detailed characters and corresponding geometrical information about the acoustic black 
hole by observing quasinormal ringdown waves in the low temperature condensed matters. 

I. INTRODUCTION 

Black holes are one of the most fascinating objects for theoretical physics, especially for General Relativity. Strik- 
ingly, they have thermodynamic properties based on the Bekenstein-Hawking entropy formula and the Hawking's 
renowned effect H, which might be a touchstone of quantum theory of gravity. Unfortunately, however, it seems to 
be impossible to detect the Hawking radiation from realistic astrophysical black holes because the Hawking tempera- 
ture of astrophysical black holes is so low that its intensity is too weak to be detected. For the purpose of confirming 
the Hawking radiation, it will be greatly helpful to perform laboratory experiments using some kinds of analogy. 

There are a lot of works about an acoustic "black hole" as the experimental field for Hawking radiation. The 
acoustic black hole has a sound horizon corresponding to an event horizon of a black hole spacetime. Unruh 
originally suggested a hydrodynamic analog of an event horizon '5| as a more accessible phenomenon which might 
shed some light on the Hawking effect. In principle, an event horizon for sound waves appears wherever there is a 
closed surface through which a fluid flows inwards at the speed of sound. In the vicinity of the surface, the flow 
is subsonic on one side and supersonic on the other side; there is a close analogy between sound propagation on a 
background hydrodynamic flow and field propagation in a curved spacetime. Thus, if such an analogous black hole is 
materialized, it will be experimentally powerful tool for confirming the Hawking radiation. There are some works in 
this direction .4, 5]. Recently, the subject has attracted much attention 6, 7, 8, 9, 10., llj . 

In order to detect the Hawking radiation in the experimental laboratory, thermal fluctuations will become main 
obstacles. Therefore we consider an acoustic black hole in the Bose-Einstein condensation (BEG) system at low 
temperature. Thanks to recent techn olog ical advances, the BEG has been observed in 1995 in a remarkable series 
of experiments on vapors of rubidium |l2j| and sodium [T^ in which the atoms were confined in magnetic traps and 
cooled down to extremely low temperature, of the order of fractions of microkelvins. (See the review paper flil about 
the phenomenon of Bose-Einstein condensation of dilute gases from a theoretical perspective.) Garay et al. jl5l ll^ 
have analyzed in detail the case of a condensate in a ring trap, and proposed a realistic scheme for adiabatically 
creating stable acoustic black/white holes. It seems to be possible to materialize an acoustic black hole in the BEG 
system. 

Astrophysical black holes will obey the uniqueness theorem for black hole spacetimes, and there will be no hair 
other than mass, angular momentum and electric charge. On the other hand, it will be possible to realize analogous 
black holes having many hairs. Since analogous black holes have geometrical correspondence, it means that many 
kinds of geometry will be realized in the analogous systems. Thus, it is necessary to identify which corresponding 
classical geometry is realized in the experiments. 

Geometries of acoustic black holes are decided by physical situations (e.g. the equations of state, boundary condi- 
tions on fluid systems and so on ). Since what directly affects sound waves or Hawking radiations is its background 
geometry or wave equation on the analogous curved spacetime, it is needed to confirm whether the required geometry 
is materialized or not by detecting some kind of a physical phenomenon which is not determined by factors except 
for the geometry. 

Without this confirmation, we can not clean away the fear of false geometries or false analogies. 

One of the ways to do it is to detect damped oscillating waves of fluid fluctuations under some perturbation. We 
make a point that the modes of damped oscillations in the linear phase are uninfluenced by given perturbation and 
directly depend only on the geometry. 

This damped oscillating waves are well known as quasinormal ringdown waves for astrophysical black holes context. 
The first order perturbed Einstein equation gives a wave equation, and quasinormal modes (QNMs) are complex 
frequency solutions of the wave equation with a purely outgoing-wave boundary condition at infinity and ingoing- 
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wave at the horizon, with vanishing incoming-wave ampHtude. They have geometrical information outside the horizon 
and, in fact, the quasinormal ringdown waves have astrophysical significance and are searched by using the data from 
the gravitational wave detectors [13 • In this paper, we analyze quasinormal ringdown waves of an acoustic black hole 
with the procedures which are frequently used in analysis about those of astrophysical ones, and using the QNMs of 
this wave, we discuss the distinguishability of the geometry. 

On the other hand, there is a subtle problem concerning the detection of QNMs in gravitational astrophysics. For 
the purpose of designing the gravitational wave detector or making the data analysis, we must predict theoretically 
amplitudes of gravitational waves which we want to detect. The amplitude depends on the energy which is extracted 
by the gravitational waves from sources such as supernovae or neutron star mergers. Thus it is important to know how 
many energies of the process transport to the gravitational waves. When it comes to quasinormal ringdown waves, 
however, it depends on the final phase of the gravitational wave emission process such as black hole formation; hence 
we need to clarify all the emission processes including nonlinear phases (e.g. the merger phase of binary neutron star 
systems) for the estimation of the amplitudes. Generally, this theoretical estimation can be done only with complicated 
numerical works which is not yet sufficiently accurate at the present time. Some data analysis has been done on the 
assumption that the energy transport is a few percent to the background black hole mass, however this assumption 
does not have much foundation |l3- If analogous black holes are realized in some experimental laboratories, we will 
obtain some suggestion about the energy extraction rate by QNMs. Therefore the study of QNMs in acoustic black 
holes is also significant in this direction. 

In this paper, we consider an analogous black hole in a Bose-Einstein condensate and predict the quasinormal 
ringdown waves in the system for the purpose of investigating their corresponding classical geometry. The BEG 
system has some remarkable properties that it is a quantum mechanical system and its dynamical evolution obeys 
the Gross-Pitaevskii (GP) equation, in which any profile of its density can be materialized freely by controlling the 
exterior potential in actual experiments. The linearized GP equation gives a wave equation which corresponds to that 
on a black hole background, giving QNMs under some appropriate conditions. In order to derive the quasinormal 
frequencies, we will use the WKB method given by Schutz and Will jl^l- Some quasinormal frequencies in the case 
where the section area of a nozzle is constant and the background density has some combinations of 'power-form' will 
be shown. This choice of the background density covers for all scale invariant forms of the density, where the power 
of the density is a unique parameter for the system. This fact will be helpful for actual experiments. 

The paper is organized as follows. In Sec. ^ we summarize the Gross-Pitaevskii equation which is the equation 
of motion of a Bose-Einstein condensed matter, and the linearized Gross-Pitaevskii equation. And then, an acoustic 
black hole which has an event horizon for sound wave is shown by considering the differential equation. The detailed 
analysis has been discussed by G. Barcelo, S. Liberati and M. Visser j23| and we summarize it in the Appendix. For 
this black hole, we calculate the QNMs by using the WKB approach in Sec. lIIII Finally, we show how to confirm the 
geometrical information and discuss the remaining problem in Sec. IIVI 



Now, in order to consider a Bose-Einstein condensate, we introduce the field operator ip which is divided into two 
parts as 



where 5* is called as the order parameter and '0 is some perturbation around the order parameter. The order parameter 
is a complex function defined as the expectation value of the field operator, and written as the following form: 



Here, the functions hq and ip denote a condensate density and a phase factor of the condensed matter, respectively. 
The order parameter is a classical field called the "wave function of the condensate" . In the following, we ignore V' 
and consider only this order parameter. 

When the number of atoms is large and the atomic interactions are sufficiently small in the vicinity of zero tem- 
perature, the Gross-Pitaevskii equation is the equation of motion of a condensed matter 



II. GROSS-PITAEVSKII EQUATION AND ACOUSTIC BLACK HOLE 



iP = 'qj + iP 



(2.1) 




(2.2) 



where m is the mass of the atoms, a is the scattering length and the atoms are trapped by an external potential 
Vext (x) . This has the form of a "nonlinear Schrodinger equation" . Substituting Eq. (|2.1ll into the above equation. 
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we obtain two equations from the real and imaginary parts, respectively: 



^ = -V(noV$), 

a$ 1, ,n 



Am? 



no 



2"o 



cxt 

m 



-no . 



Here we can introduce a new function $ = , and the velocity of a condensed matter vg as 

vo = V$, 

and above two equations 1)2. 3() and 1)2. 4|l seem to be the continuity and the Bernoulli equations, respectively. 
Next, we derive the equation for the linear perturbation of the Gross-Pitaevskii equations (|2.3(l and (|2.4(l 

no = no(l + C) , 
$ = $ + ?7 . 



(2.3) 
(2.4) 

(2.5) 



(2.6) 



Here we consider no and $ as a background, and C and rj represent small perturbations around them. For simplicity, 
we assume that the background is stationary and independent of time. Then, the background values no and are 
obtained by solving the following time-independent Gross-Pitaevskii equation: 
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V^no (Vno)2 
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2 ^ ' Am? \ no 
On the other hand, the linear perturbation equations are given by 

^ -X/t A-\7'^ 

drj _ f? 1 ^, ^ ^, An af? 



V(no V*) 

Airah? 




0, 
0. 



dt ^ ^ Am"^ TiQ 



vo • VC + V^?y + Vlnno • V77 = 0, 
V(noVC)+ ''17' noC = 0. 



(2.7) 
(2.8) 

(2.9) 
(2.10) 



Now we assume that the variation of the density no can be neglected for the healing length (= 1/ ^/iirano), i.e., the 
third term on the L.H.S. of Eq. (|2.1()|l is neglected. This approximation is considered as a low frequency approximation. 
And then, we obtain the simple equation: 

Airaf? ^ , , 
^noC - 0. (2.11) 



drj 

-r- + vo • Vt? 

at 



Eliminating the density perturbation C, from Eqs. H2.9|l and H2.11|l . we obtain a second order differential equation for 
the phase (velocity) perturbation -q 
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^ + vo • V U 



dt 



— V(noV?7) = 0, 
no 



(2.12) 



where we have defined as the sound velocity: 



Cs — — \/4 TT a no . 
m 



(2.13) 



It should be noted that the sound velocity Cg is a local quantity and also that if this low frequency approximation is 
not taken into account, Ea. H2.9|l and (|2.1U|) can not be reduced to a single second order differential equation for the 
condensate phase perturbation. We will give some discussion about this approximation in Sec. IIVI 

In the following, for simplicity, we consider the case where the background velocity has only x-component, i.e., 
{vq} = {vq, 0, 0} where the superscript i denotes {x, y, z}. From Eq. H2.12|l . we can read off the distance between 
two nearby points in the corresponding curved spacetime as 

(~(Cs ~ VQjdt — 2vQdtax + dx" 



ds' = 



auQ 

Cs 
Cs 



-{cl 



dt 



rdx 



dif + dz^) 
^2 



^dx^ + dy'^ + dz' 



- 1- 



2 J 2 
2„ dr 



-dx^ + dy'^ + dz'^ 



(2.14) 
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This detailed derivation is shown in the Appendix. Here, we have introduced the new time coordinate r defined by 

dT^dt+ „ dx. (2.15) 

When the superfiuid velocity changes from subsonic to supersonic velocity, there is a sound horizon for the direction 
of the X-axis. We shall call this system with the horizon an acoustic black hole, because this has an event horizon for 
sound waves. 



III. QUASINORMAL MODES 

In this section, we consider the background configuration of a condensed matter which has a sound horizon and 
then derive the quasinormal frequencies in this background. In order to obtain a sound horizon, we consider the one- 
dimensional background flow of a condensed matter with the velocity vq in the x-direction. The background density 
uq is derived by using the continuity equation (|2.7|l . Then, we can calculate the external potential by Eq. (|2.8() . In 
practice, if the external potential is chosen as the above value, it is possible to construct an acoustic black hole. We 
note that the background discussed in the following has a constant section area of a nozzle. This is a single black hole 
case. 

Using Eq. I|A1|) . the linear perturbation equation of the system is written as 



no 



no vo 



+ (^no(^l-^)a. 



= 0. 



(3.1) 



For the above equation, using the background density no, we rewrite the sound velocity Cs as 



(3.2) 



and the background velocity vo can be expressed as 



V 

vo = — . 
no 



(3.3) 



Here Eq. (|3.3|) is obtained from the continuity equation (|2.3|) and the assumption that the section area of a nozzle is 
constant. Then we rewrite the above perturbation equation into the following form: 



2V , 



no 



1 



no 



-at dt d.,^v[ d,— dt + c'd, „o n _ a 



V' 1 



C2 n3 



77 = 0. 



(3.4) 



We assume that the background is static, so it is convenient to consider the Fourier coefficients of the phase (velocity) 
perturbation ry 



X^{x) = / dtr]{t, x)e' 



(3.5) 



In order to simplify the perturbation equation, i.e., to remove the first differentiation with respect to i, we introduce 
the following auxiliary functions, F^^ and G^: 



(3.6) 



where the function F^^ has been chosen as 

F^ix) 



dx ■ 



ujV 



(3.7) 



Then, the perturbation equation takes the following form: 



c -(- f. Ox 



(3.8) 



where we have defined (3 = V/G and a = lj /G , respectively. 
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Next, in order to obtain the Schrodinger-type equation, we introduce the new coordinate as 

da;, 1 



5/2 

From the above coordinate transformation, the perturbation equation becomes 



dx /(a;)' 

m = ^i^. (3.9) 



1 {d^no)Jj4-j3^ 2 

^x, "t" f, 7/2 '^x, -r o 



G^=0. (3.10) 



Furthermore, in order to remove the first differentiation with respect to a;,, we replace the function G^j as 

G^{x) = n-^'^H^{x). (3.11) 
Then we obtain the Schrodinger-type equation for the function H^j as 

" 4 < 16 ^ • ^ 

In order to calculate the quasinormal frequencies, it is necessary to choose the appropriate potential V. For example, 
we consider the case that the background density uq has the following form: 

no = Nx" , (3.13) 

where a is a constant. From Eq. 1)3. 3|l . the background velocity becomes 

vo = ^x- . (3.14) 

It is noted that the different choice of a derives the different relative velocity profile in the fluid flow. In this case, the 
horizon arises at 

V) . p^i^) 

because of the sound velocity Cg = C\/]Vx°/^. 

We flnd that the range of parameter a should be —2/5 < a < 2 for the condition in which the potential V vanishes 
at infinity. Furthermore, in order to obtain the potential which goes to zero at the horizon of an acoustic black hole, 
it is necessary to restrict a to < a < 2. This condition for a also ensures the short-range character of the potential 
V which has been mentioned in 

V(x,)(ix, < CO. (3.16) 

Thus, in the following, we consider the case of the parameter range < a < 2. To obtain the above background 
density no, we may choose the external potential Vcxt of the background Gross-Pitaevskii equation as one calculated 
by using Eq. (|2.8() . In the following discussion, we choose (3—1 and iV = 1. For example, in the case of a = 5/3, 
substituting no given by Eq. H3.13|l into the perturbation equation l|3.12(l . we obtain the simple potential V and the 
coordinate transformation / as 

5 (a;5 - 1)(2;5 -I- 19) 
48 ^^3175 ' 

a;^ — 1 

f = ^- (3-17) 
The potential V shown in Fig. ^ has one peak and goes to zero at a;, — > ±oo in the a;, coordinate. 



FIG. 1: The function V{x) in the case of no = x^^^. V{x) = at a:: = 1 (horizon) and the Hmit x — > oo. 
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5R(a) 


9(a) 





1.075 


-0.796 


1 


1.685 


-1.523 


2 


2.132 


-2.006 


3 


2.501 


-2.394 



TABLE I: Some quasinormal frequencies for a = 5/3. 



In order to obtain the quasinormal frequencies of the background with the above potential, we use the usual WKB 
approach. (See Schutz and Will [l^ in the case of a black hole. The higher order WKB approach has been developed 
in Refs. [2lll2^l2^ .) The quasinormal frequencies are obtained as 

a^ = Vo-i(n+^Y-2Vo)^/\ (3.18) 

where n is the non-negative integer for ^{a) > 0, n is the negative integer for ^(a) < 0, and the prime ( ) denotes 
the differentiation with respect to . Here n indices the n-th QNM, and the subscript on the potential denotes the 
value at 

a;o = 1.152 where Vg = . (3.19) 

The differentiation of the potential in the above calculation are obtained directly. The quasinormal frequencies are 
shown in Table I. 

Now, we have a interest in the power a dependence of the quasinormal frequencies. The quasinormal frequencies in 
the cases of a = 1/3, 2/3, 1, 4/3 and 5/3 are summarized in Fig. [3 The figure shows the least damped mode in each 
value of a differs each other, and thus, we can distinguish the power a by the quasinormal frequencies and check the 
realization of the corresponding geometry in the analogous black hole system. 

In Fig. 13 the QNMs in the case of no = (6 a; -I- (1 — b)x^y^^ are shown, where the parameter b takes the value 
of 0, 1/3, 2/3 and 1. In this case, the sonic horizon exists at the same position {x = 1) as in the previous case of 
no = x". We find that as 6 1, the QNMs change into those of no = x^^^ continuously. 



IV. DISCUSSION 



In this paper, we have discussed the QNMs of an acoustic black hole in the context of the Bose-Einstein condensate. 
This acoustic black hole is prepared by choosing the external potential, and a single black hole case has been discussed. 
Using the WKB method [l9| , we have derived some quasinormal frequencies in the case where the background density 
has a power form or some combination of power forms. 

The QNMs in this system can be considered as the ringdown waves emitted from a perturbed acoustic black hole as 
in the usual astrophysics context. We can obtain the detailed characters and corresponding geometrical information 
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FIG. 2: The quasinormal frequencies in the case of no = x" for various o. The symbols o, +, □, o and O show the case of 
a — 1/3, 2/3, 1, 4/3 and 5/3, respectively. 
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FIG. 3: The quasinormal frequencies in the case of no = (6 a; + (1 — 6) 
& = 0, 1/3, 2/3 and 1, respectively. 
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The symbols Q, 



and o show the case of 



about the acoustic black hole from its detection as follows. It is our purpose mentioned in Sec. ^ From the quasinormal 
frequencies shown in Fig.|21and Fig.|21 we can identify the density of the BEC which determine the potential function 
(|3.12l) . Thus, they are useful to check whether the system is realized to our expectations in an actual experiment and 
to know the potential function (|3.12(l and the corresponding geometrical information such as Eq. I|A4|I . Not only the 
least damped mode but the higher damped modes are useful to identify the corresponding classical geometry. 

It is noted that |3fi((T)/3((T)| ~ 1. This means that the ringdown waves from this system have a short decaying 
time. In astrophysical black holes, the situation is similar: their ringdown waves have a short decaying time. One 
exception occurs in the case of extremal black holes whose Hawking temperature equals zero. Thus, at the present 
moment, we do not have a strong interest in such black holes. 

We also note that "Schwarzschild like black holes" can not be obtained in the fluid configuration discussed in this 
paper. However the Schwarzschild black hole is one of the most impotant objects in astrophysical black holes context. 
If the acoustic black hole which has the same geometry as the Schwarzschild black hole is realized in some experiment, 
it is quite meaningful. Why we could not obtain the geometry is because the sound velocity Cg and the background 
velocity vq are related to the density no as E g. (13.311 . In order to obtain a Schwarzschild black hole in the BEC, 
we have to consider other fluid conflgurations |2J|. And then, the quasinormal frequencies are solvable by using the 
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Leaver's method |25j in which we may solve the three-term recurrence relation. 

In our study, we have considered the low frequency approximation for the linearized Gross-Pitaevskii equation. 
This is relevant in the region where the wavelength of sounds is much longer than the healing length of the system. 
The quasinormal ringdown waves obtained in this paper satisfy the condition and therefore, our investigation is 
appropriate. 

The QNMs obtained in this paper depend on the outer geometry from the sound horizon. The higher damped modes 
which are not obtained in this paper will depend on the geometry in the vicinity of the horizon, where, however, the 
wavelength becomes shorter and shorter due to the blue shift effect associated with the horizon. When we intend to 
obtain the higher damped modes and, correspondinly, the information about horizon, we cannot use the low frequency 
approximation and have to investigate the full linearized Gross-Pitaevskii equation. 

In this case, we cannot obtain a single second order differential equation for the condensate phase perturbation, and 
sound waves will break the causality because th eir prop agation speed may become infinity. There are many papers 
about this "broken" Lorentz invariance problem |26l Wa |29| . 

Recently, Konoplya studied the influence of the back reaction of the Hawking radiation upon QNMs in the case of 
(2-1-1) dimensional black hole 'so]. Such effects might be more important in this BEG system because the influence 
will give some information about the Hawking radiation at low temperature. 

We leave these issues for future investigations. 
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APPENDIX A: ANALOGUE GRAVITY IN THE BOSE-EINSTEIN CONDENSATION 



In this appendix, we summarize the derivation of the curved metric (|2.14|l from the wave equation H2.12|l . More 
detailed analysis has been discussed in poj . 

Using the background continuity equation H2.7|) . we rewrite the wave equation (|2.12|) for the linear perturbation of 
a velocity potential as 



1 fd 



no \dt 



+ (V • vo) 



no f d 



1 \dt 



Vo • V ?7 



no 



-V(noV77) = 0. 



(Al) 



Now, let us compare the above equation with the scalar wave equation in a general curved spacetime: 



I d_ 



— d 



(A2) 



where g^'^ denotes a metric and g is its determinant. We choose the coordinate as {x^} — {t,x,y,z}, and if we set 
the metric components as 



no 



\ 



-99'^ 



no Wo 



V 



(A3) 



where the superscript i denotes {x, y, z}, then the wave equation ljA2p is equivalent to (jAl|) . Using the value of 
determinant y^—g = a^no/cg, the inverse matrix is calculated as 



ano 




(A4) 
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Thus we can interpret that the perturbed system characterized by Eg. ljAll) corresponds to the wave propagation on 
the curved spacetime background with the metric (|A4() . 
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